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We use the symmetries of monolayer graphene to write a set of constraints that must be satisfied 
by any electron-phonon interaction hamiltonian. The explicit solution as a series expansion in the 
momenta gives the most general, model-independent couplings between electrons and long wave- 
length acoustic and optical phonons. As an application, the possibility of describing elastic strains 
in terms of effective electromagnetic fields is considered in detail, with an emphasis on group theory 
conditions and the role of time reversal symmetry. 
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I. INTRODUCTION 

Symmetry plays an obiquitous role in solid state 
physics, where invariance of the system under symme- 
try transformations imposes strong constraints on the 
form of hamiltonians and other physical observablesji*^ 
In computations of phonon spectra, the choice of conve- 
nient symmetry-adapted coordinates simplifies the diago- 
nalization of the dynamical matrix.'^ Instead of a 3iV x 3 
matrix, where N is the number of atoms in the unit cell, 
we end up with a collection of smaller matrices, one per 
inequivalent IR (irreducible representation), with a di- 
mension equal to the multiplicity of that representation 
in the 'mechanical' representation. Similar simplifica- 
tions take place in the computation of electronic bands. 

An important caveat is that the symmetry which is 
effective in simplifying the problem is not the whole 
space group of the crystal, but the little group^ that pre- 
serves the wavevector of the electrons or phonons J' Thus, 
symmetry-based methods are really powerful around 
points and lines of high symmetry in the Brillouin 
zone. Since electron-phonon interactions involve generic 
wavevectors with trivial little groups, it would seem 
that symmetry considerations can not be very helpful 
in that case. Indeed, many computations of electron- 
phonon interactions are based on the use of specific mod- 
els (electron-ion interactions,'^'^ modulated hoppings^i^i^ 
phonon modulated electron-electron interactions^^ etc), 
and little or no use is made of the symmetries of the 
system. 

On the other hand, long-wavelength phonons play a 
crucial role in many physical processes of interest, such 
as transport phenomena, where electron-phonon inter- 
actions are essentiali^i^ The elements of the dynamical 
matrix for long wavelength phonons (acoustic and opti- 
cal) can be expanded as a power series in the phonon 
wavevector q, and this expansion is strongly constrained 
by the little group of the T point, which of course coin- 
cides with the space group of the crystal. Unfortunately, 
low energy electrons lie near the Fermi surface, where 
typical points lack any symmetry. 

Grapheneiiii^ is a remarkable exception in this re- 
gard. Instead of an ordinary Fermi surface, undopped 



monolayer graphene has two isolated, high symmetry 
Fermi points. Thus, as long as we restrict ourselves to 
long wavelength phonons and low energy electrons, which 
are necesarily close to the Fermi points, the interaction 
hamiltonian will be strongly constrained by symmetry. 
In fact, in Section 4 we will be able to write the most 
general interaction hamiltonian as an explicit series ex- 
pansion around the F point (for phonon wavevectors) and 
the Fermi points (for electron wavevectors). 



One possible application of the hamiltonian obtained 
here is to many-body computations,— i^^'^"^ where analytic 
expressions are often preferred over numerically defined 
functions. In the usual approach a concrete model, such 
as one based on phonon-modulated nearest neighbors 
hoppings, is used to compute an interaction hamiltonian 
which later is expanded to the required order in the mo- 
menta i^°i^^'^^ Alternatively, one can take the interaction 
hamiltonian given in Section 4 as a power expansion, and 
keep terms up to the order desired. Obviously, less effort 
is involved by this method, and one is sure that all pos- 
sible couplings have been included, whereas a particular 
model may miss some of then. The obvious drawback is 
that symmetry arguments by themselves do not give a 
clue as to the values of the parameters. Another, very 
different application of the electron-phonon hamiltonian 
is presented in Section 5, where we investigate a descrip- 
tion of elastic strains in terms of effective electromagnetic 
fields^ 



This paper is organized as follows. In Section 2 we 
review the symmetries of monolayer graphene and the 
transformation properties of electrons and phonons. Sec- 
tion 3 presents a study of the dynamical matrix includ- 
ing the form of the long wavelength normal modes. The 
phonon-electron interactions are computed in Section 4, 
and the results are used in Section 5 to study the possi- 
bility of interpreting trilinear electron-phonon couplings 
in terms of effective scalar and vector potentials. The 
conclusions are presented in Section 6. 



II. SYMMETRIES 

We begin this section with a description of electrons 
and phonons in monolayer graphene, including the choice 
of the symmetry-adapted modes to be used in the rest of 
this paper. We then proceed to study their transforma- 
tion properties. 



A. Electrons and phonons 




Monolayer graphene consists of a honeycomb lattice. 
The corresponding hexagonal Bravais lattice is gener- 
ated by ti = (a\/3, 0) and ^2 = (a\/3/2, 3a/2), where 
a is the distance between nearest neighbours. (See 
Fig. [1]). The unit cell contains two atoms, with positions 
given by ri =ii/3 -1-^2/3 and 7^2 — 2ri. The recipro- 
cal lattice is generated by gi — (27r/V3a, — 27r/3a) and 
52 = (0,47r/3a), with U ■ gj = 2T:dij. Corresponding to 
the two atoms in the unit cell, we may define two Bloch 
wave functions 



i) , z = l,2 (1) 



where the sum runs over all the points in the direct lat- 
tice, i.e., t = niti -\-n2t2 and $(r) is a real 7r-type atomic 
orbital. As is well known, a simple tight-binding com- 
putationiiii^ yields a spectrum with two Fermi points 
located at Ki = — 2gi/3 — 32/3 and K2 = —Ki. Near 
the two Fermi points, the hamiltonian can be linearized 
and one finds 



(2) 



where t is the hopping integral, and H{—Ki + fc) = 
H*{Ki — k). Thus, the low energy electronic excitations 
behave like massless Dirac fermion a^^'^° with relativistic 
spectrum E = ±Vi?|fc|. Henceforth we assume our units 
are such that vp = 3at/2 = 1. 

In this paper we want to consider the interaction be- 
tween electrons near the Fermi points {Ki, —Ki} and 
phonons Q{q) for small values of q. Due to invari- 
ance under reflection by the horizontal plane cr^, only 
in-plane modes can couple linearly to electrons. Thus, 
for most of this paper we consider only in-plane modes, 
which are specified by giving the horizontal displace- 
ments u{t ) = (cci, yi, X2, 1/2) of the two atoms in the unit 
cell with origin at the lattice point t. See however Sec- 
tion 5, where quadratic couplings of out-of-plane acoustic 
modes are considered in the long wavelength limit. 

We can choose a basis of symmetry-adapted modes 
with particularly simple transformation properties by 
noting that, under a point group transformation, a mode 
undergoes two changes: The vectors describing atomic 
displacements are rotated or reflected according to the 
vector representation, and there is a permutation among 
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FIG. 1: (Color online) Top: Direct lattice and unit cell for 
monolayer graphene. Bottom: First Brillouin zone and Fermi 
points. The vectors K'^,K'( {K '2,1^2) are equivalent to K\ 



equivalent atoms. This last effect may give rise to an- 
noying phases in the transformation matrices, which can 
be avoided by incorporating appropriate phases into the 
definitions of the modes. The use of circular rather that 
linearly polarized modes simplifies the effect of rotations. 

Our choice of in-plane modes is given in TableU In the 
long wavelength limit qa « 1 , the phases vanish a {q) — > 
1 and the unprimed (primed) modes become circularly 
polarized acoustic (optical) normal modes respectively 
(see Section 3 for details). 
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TABLE I: In-plane symmetry-adapted modes for monolayer 
graphene. Here 6i(g) — ^e^'''^^ , where {ri} are the equilibrium 
positions of the two atoms in the unit cell. 

This basis is used to expand general in-plane displace- 
ments u(t) = (a;i, j/i, a;2, 2/2) as follows 



q \i=l 1=1 / 



(3) 
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B. Transformation properties 

The space group of monolayer graphene is symmor- 
phic, with point group Dq^ (see refi^ for conventions and 
notation). Dg/j is generated by {C^, u^i, cr/j}, where Cq 
is a counterclockwise rotation by tt/B around the z axis, 
(T„i is a reflection by the xz plane and ah is a reflection by 
the horizontal xy plane that contains the graphene layer. 
However, given the planar structure of the system and 
as long as we are restricted to in-plane modes, we will 
consider instead the normal subgroup Cqv, generated by 
{CQ,avi}. The little point group of Ki is — C^v, 
generated by {C^,avi}- The star of Ki consists of the 
two Fermi points {Ki, —Ki}. Since the binary rotation 
C2 around the z axis transforms Ki into --fCi, it is con- 
venient to use {Cg", a^i, C2} as a (redundant) set of gen- 
erators for Cqv 

The action of the generators of Cqv on the modes in 
Table I is given in Appendix A. As a result of a sym- 
metry operation g, a displacement with wavevector q is 
transformed into a new displacement with wavevector 
gq and coordinates {gQi,qQ'i}- Defining Q = {Qi,Q2} 
and Q' = {Qi, Q'2}j the transformation properties can be 
written 

gQi9q)^Tig)Qiq) , gQ'igq) ^T'{g)Q'{q) (4) 
where T'(C2) — — 7(6*2) = 1 is the unit matrix, and 

1 



(5) 



with aj = exp(27ri/3). 

Another important symmetry is time reversal, which 
acts by complex conjugation. According to Table I, 



T : ul{q) ^ U2{-q) , u^iq) ^ ui{-q) 



(6) 



with identical action on u'^. Then, reality of a general 
displacement ^ implies 



'{q)^<y.Q'{-q) (7) 



We now turn our attention to the electrons. Since 
low energy electronic excitations are located near the 
two Fermi points, it is convenient to change variables, 
K — ±Ki + k. Correspondingly, we define Bloch func- 
tions $i± (fc) in the neighbourhoods of the Fermi points 



$,±(fc) = $,(±i?i + fc) 



(8) 



The action of the generators of Cgi, on the wavefunc- 
tions is given in Appendix A. As long as we do not con- 
sider the spin, time reversal acts simply by complex con- 
jugation. According to ([T]), $*(i?) = (f>j(-i?). Thus 
time reversal exchanges the two Fermi points and we have 



r : $,+ (fc)-><i>,^(fc) = $,_(-fc) 

<i>,_(fc)-><f:_(fc) = $,+(-fc) 



As both time reversal T and the binary rotation C2 
around the z axis reverse the sign of the wavevector, 
the combined symmetry TC2 preserves the Fermi points. 
Note however, that the two atoms are exchanged 

TC2 : $i+(fc) ^ $2+(fc) , $2-(fc) ^ $i-(fc) (10) 

This combined symmetry plays an important role in the 
computation of electron-phonon interactions in Section 4 
and is behind the topological stability of Fermi points in 
mono and multi-layer graphene 1^ 



III. DYNAMICAL MATRIX 

The modes in Table I have simple transformation prop- 
erties, but they do not correspond to normal modes with 
well defined frequencies. In general, these must be ob- 
tained by numerical diagonalization of the dynamical ma- 
trixi^i^i2^ However, in the long wavelength regime one 
can easily obtain useful analytical results by exploiting 
the symmetries of the system. 

The symmetries constrain the form of the dynamical 
matrix V{q) defined by 



(11) 



where Q = {Q,Q'} stands for the four symmetry- 
adapted coordinates in (l3|), m is the mass of the carbon 
atom and the sum is over the first Brillouin zone. The 
hermitian matrix V can be decomposed into 2x2 blocks 



nq) = ( ^^t D' 



(12) 



associated to the unprimed and primed modes. Then, 
invariance under a spatial symmetry ([4]) gives 



T{g)V{q)THg) = Vigq) 



with 



T(q) ( ^(-9) 



whereas invariance under time reversal ([7]) implies 

CJ:r. 



TV*{q)T^V{-q) , T = 



CT. 



(13) 



(14) 



(15) 



Neglecting initially the acoustic-optical mixing due to 
M, we can obtain an approximation to acoustic frequen- 
cies and modes by solving a 2 x 2 eigenvalue problem 



D{q)u{q) ^ uj\q)u{q) 



(16) 



(9) 



where u{q) is a linear combination of the two unprimed 
modes given in table I. As shown in Appendix B, the 
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symmetry constrains the form of the hermitian matrix D 
to 



D{q) 



aiq) b{q) 
b*{q) a{q) 



(17) 



where a{q) and b{q) are even functions of q that satisfy 

(18) 



a{ujq) — a{q) , b{Loq) = u}b{q) 
a{q*) = a{q) , b{q*) = b*{q) 



The eigenvalue problem is solved by 

ujl{q)=a{q)±\b{q)\ 
with normal modes given by 

1 . , , biq) 



u±{q) 



±U2iq)) 



(19) 



(20) 



In the absence of acoustic-optical interactions, this is 
the most general acoustic spectrum compatible with the 
symmetries of the system. For small q this solution ad- 
mits a simple geometric interpretation. Expanding a{q) 
and b{q) 



aiq) = a2\q\'+ai\q\^ + 0{q'') 

b{q) = 629*' + 64|gp9*' + M' + 0(q6) 



(21) 



where all the constants are real, and keeping only the 
quadratic terms yields 



^Lriq) = («2 ± &2)|gr 



(22) 



where the + (— ) sign goes with Luj^ (w^). The normal 
modes are 



UL{q) 



1 



V2\c 



UT[q) = - 



x/2| 



91 



q*ui{q) +qu2{q) 



q*ui{q) - qu2{q) 



(23) 



These represent longitudinal (LA) and transverse acous- 
tic (TA) modes respectively. To see this, note that for 
qa « 1, we have 



y-hiq) 




UT{q) 



1 

71 



/ — sin 6 

cos 6 



— sm f 
cos 6 



(24) 



where q 



\qW\ i.e. 



= ia.n~^ {qy / qx) ■ In the long 
wavelength limit, the normal modes are completely inde- 
pendent of the coefficients in the dynamical matrix. 

We can easily go beyond the quadratic approxima- 
tion in q valid in the long wavelength limit by keeping 
more terms in the expansions pTjl . Keeping just the first 
anisotropic contributions gives 

^Iriq) = (a2±&2)k|' + |<z|^a4±&4±&4Cos60) + O(q6) 

(25) 



One can use perturbation theory to compute the leading 
corrections ul,t — ^ ul,t + Sul,t to the normal modes. 
Taking as perturbation the fourth order corrections to D 



SD ^ 



ai\q\ 



bi\q\'q*^ + b4q^ 



b,\q\\^ + biq*^ ai\q\' 



(26) 



yields 



Sul 
Sut 



Urp SD ul 





UJ 



\ SDuT 

2-"i 



2h 



2b2 



Igp sin 66' ut 



P sm66uL 



(27) 



Thus, beyond the long wavelength approximation the 
normal modes are not purely transverse or longitudinal, 
and (p7|) shows that this longitudinal-transverse (L-T) 
mixing is an 0{q^) anisotropic effect. 

The solution to the eigenvalue problem for the opti- 
cal modes is very similar. The elements of the matrix 
D'{q) satisfy the same constraints (|T8l) . but the diagonal 
elements need not vanish for q = 0. Keeping just the 
quadratic contributions 



D'{q) = 



4191' 



b'2q' 



: + a'2\q\'- 



(28) 



we find that the optical normal modes and ui^ are 
given by (|23p in terms of the primed modes u'^ , with nor- 
mal frequencies 



^ZriQ) =^0 + (4 ± b'2)\qf 



(29) 



Similarly, L-T mixing for optical modes is still given 
by ((27|) if one uses primed constants. 

We can also study acoustic-optical (A-0) mixing, 
which is induced by the matrix M. For example, keeping 
only the leading contributions to M (see the Appendix B 
for details.) 



M{q) ~ ei 



q 
-q* 



(30) 



and using perturbation theory yields the following cor- 
rections to the LA and TA modes 



Sui 



Sut 



lei 



lei 



gl (sin 36 u'j^ + cos 36* uip) 
(/I (cos 39 u'l^ — sin 30 u'rp) 



(31) 



Thus A-0 mixing is an 0{q/uj'^) anisotropic effect. Ex- 
changing primed and unprimed modes in (j3ip gives the 
corrections to the optical modes. This mixing induces a 
common shift in the frequencies of the acoustic modes 



(32) 



while the opposite shift is induced in the optical frequen- 
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IV. ELECTRON-PHONON COUPLINGS 

In this section we use a similar strategy to obtain the 
most general interaction hamiltonian compatible with 
the symmetries of the system. We consider an electron- 
phonon interaction hamiltonian of the form 

H.-ph = Q)4+(fc + 9/2)c,+ (fc - q/2) 

ij,k,q 

+ ^ iJ_,,(fc,g,Q)4_(fc + g/2)c,_(fc-g/2) (33) 

ij,k,q 

where Ci±{k) annihilates an electron in the state 4'i±(fc) 
given by ([8]) and Q = {Q, Q'} stands for the four 
symmetry-adapted coordinates in ([3]). The somewhat 
unusual parametrization of initial and final electron mo- 
menta simplifies the form of the constraints. 
The matrix element is given by 

iJ,,(fc,g,Q) = ($,_,(fc-f|)|F,,|ci>^.+ (fc- |)) (34) 

where Vei is the electron-ion potential. Since the two 
Fermi points are related both by time-reversal and by 
C2, either symmetry can be used to obtain the form of 
the hamiltonian near K2 = —Ki, giving 

H^{k,q,Q) ^ H*{-k,-q,Q) (35) 

where 

iJ_,,(fc,(7,Q) = ($,_(fc + (7/2)|14»|*j-(fc-g/2)) (36) 

Thus, we can restrict ourselves to the hamiltonian in the 
vicinity of Ki . 

On the other hand, the combined symmetry C2T (jlOp 
leaves the two Fermi points invariant and can be used to 
impose a constraint on H{k, q, C2Q) 

H(k,q,C2Q)^a,H*{k,q,Q)<^x (37) 
I 



The functions entering the hamiltonian must satisfy 
some simple conditions. In particular the functions /, 
/', §2, and 32 must satisfy identical constraints, which 



which is solved by 

n(hnC^\-( Oiik,q,Q) P{k,q,Q) \ /oo^ 
^ U*(^''?'^2Q) a*{k,q,C2Q) J ^'^^> 

where a and /3 are arbitrary complex functions. Note 
that hermiticity of the hamiltonian does not imply re- 
ality for a. The point is that the matrix if(fc,g, C2Q) 
connects initial and final states with different momenta, 
and instead we have 

iJ,,(fc,g,Q) = (<i>,;+(fc+|)|Kd$,+(A;-|)) 
- (<f,+ (fc - |)|K,:|<i>.+ (A; + |))* = H;,ik, -q, Q) (39) 

which implies 

a{k,-q,Q) = a*{k,q,Cl) 

/3(fc,-9,Q) = f3{k,q,C2Q) (40) 
Finally, we must impose invariance under the little 
point group C^y 

T\g)H{k, q, Q)T(.g) = H{gk, gq, gQ) (41) 

where T{g) is given by ([5]). The hamiltonian p8|l to- 
gether with the constraints (|40l4ip contains the most gen- 
eral interactions beween electrons and in-plane phonons 
compatible space group and time reversal symmetries. 
These constraints are analyzed in detail in Appendix C, 
where the general solution for the interaction hamilto- 
nian is found in the form 

Hik, g, Q) = Hik, g, Q) + H\k, q, Q') (42) 

with 



(43) 



we write only for / 

f{u:k,u:q)^uf{k,q) , f{k*,q*)^f{k,qY (44) 



H{k,q,Q) 



f{k, q)Qi{q) - f{k, ~qyQ2{q) 9l{k, q)Qi{q) + 52(fc, q)Q2{q) 
52(fc, qYQliq) + 9l{k, q)*Q2{q) -f{k, -q)Qi{q) + .f{k, q)*Q2{q) 



H'{k,q,Q') = 



f'{k, q)Q[{q) - /'(fc, -g)*Qi(g) g[{k, q)Q[{q) + g'2{k, q)Q'2{q) 

-g'2{k,qrQ[{q) - g[{k,qrQ'2{q) f'{k,-q)Q[{q) - nk,qYQ'^{q) 

I 
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whereas the conditions on gi and g'l are 

gi{ujk,ujq) = gi{k,q) , gi{k* ,q*) = gi{k,q)* (45) 
Besides, we must have 

9i{k,-q) ^ -gi{k,q) , g[{k,~q) = g[{k,q) (46) 

Since lo^ = 1, one immediately sees that any combina- 
tion of monomials of the form 



(47) 



with real coefficients and n~n' + m- 



= 1 mod 3 will 
satisfy (j44p . Imposing instead n — n' + m—m' = mod 3 
we will get solutions to (US)) . Finally, the conditions (HS)) 
are satisfied by taking m + m' odd (even) for gi (y^). 

Eq. (|43|) together with the comments around eq. ((47)) 
can be used to expand the interaction hamiltonian to 
any order in q and fc, giving all possible in-plane phonon- 
electron couplings compatible with the symmetries of 
monolayer graphene, and are the main result in this pa- 
per. For the rest of this section we will study the long 
wavelength limit and some leading corrections. 

Consider first the coupling to acoustic phonons. Tak- 
ing into account that the interaction should vanish for 
g — > 0, we can immediately write the leading contribu- 
tions to the functions /, gi , g2 

f{k, q) « V2aiq , gi{k, g) « , g2{k, q) « V2f3iq (48) 

where ai and /3i are real constants and the has been 
introduced for later convenience. This gives 



H{k,q,Q) 



(49) 



q*Q2{q)) PiqQ2{q) 
Piq*Ql{q) a,{qQ,{q)+q*Q2{q)) 



For some physical applications it may be more conve- 
nient to express the hamiltonian as a function of normal 
rather than symmetry- adapted coordinates. As we know 
from the previous section, the longitudinal and transverse 
modes in (1231) are a good approximation to the actual 
normal modes in the long- wavelength limit. Using 



/I = ^— (Ql-zQt) , Qi^^-iQL+iQr) 



V2\c 



(50) 



the interaction hamiltonian takes the following form 

Hik,q,Q)^ (51) 
2aiQL l3ie^'\QL+iQT) 



i\q 



(iie-''%QL-iQT) 



where — tan~^(gj,/(7a;) is the phase of g= l^je*^. This 
shows that TA phonons do not couple diagonally to elec- 
trons in the long- wavelength limit, which is parametrized 
by the two real coulplings ai and This form of 

the hamiltonian can be used to write p3p in terms 



of phonon creation and annihilation operators;^ with 



QL^riq) = (2a;L,T)-'/'(aL,T(g) 



i-q))- 



Note that even if one is primarily interested in acoustic 
phonons, because of the acoustic-optical mixing discussed 
in the last section, one may have to consider the couplings 
of electrons to optical modes as well. Since the couplings 
to optical modes do not have to vanish as q ^ 0, the 
leading contributions to the hamiltonian are independent 
of q. These are found to be 



H{k,q,Q')^zV2p[ i 



Q[ 



?2 

-ie 



[Q'l - ^Q't) 




(52) 



where f3[ is a real constant. 

We can easily include higher powers of the momenta. 
For instance, for acoustic phonons the contribution of or- 
der 0{q^) is parametrized by a single real constant a2 and 
is purely diagonal and proportional to the Pauli matrix 



AH{k, q, Q) = a2\q\^{QL sin 361 + Qt cosW)a, 



(53) 



This is qualitatively different from the leading contribu- 
tion ([5T|) due to the anisotropy in the strength of the 
electron-phonon couplings and the diagonal contribution 
of TA phonons. The contribution of order 0{q) for opti- 
cal phonons is proportional to the unit matrix 



AH'{k,q,Q')^ia'2\q\Q'L^ 
and there is also an 0{k) term given by 



(54) 



Ai7'(A:, 9, Q') = aal/cKQi sin((^ - B) + Q'j, cos(<^ - 0))1 

(55) 

where f = tan ^{ky/k^) ■ The contributions of order 
0{qk) are more complicated and depend on five real con- 
stants. Note that as we move away from the F-point the 
LA and TA modes given by (|^5|) are no longer a good 
approximation to the normal modes of the system, and 
one has to take into account the L-T and A-0 mixings 
considered in the previous section. 

Electron-phonon interactions can be rewritten in terms 
of electrons and holes. The eigenstates for the electron 
hamiltonian ^ are given by 



= J^fe-"^/^^i{k) ± e''^/2$2(fc) 
V 2 ^ 



(56) 



and this can be used to obtain the matrix elements be- 
tween electron eigenstates, with and (_ff_| 

and H [_) corresponding to intraband (intcrband ) tran- 
sitions. Of course, if we decide to include higher pow- 
ers of k in the electron-phonon hamiltonian, for con- 
sistency we must also go beyond the low-energy Dirac 
hamiltonian The corrections to the electronic hamil- 
tonian are fixed by symmetry. Note, in particular, that 
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the electronic hamiltonian can be considered as the Q- 
independent part of Then (|38p reduces to 



H{k) 



(57) 



where hermiticity (|40p imphes that is reaL The httle 
group constraints (jCip simphfy to 

Cg^ : ao(a;fc) = ao(fc) , ^o(^^fc) = ^^*/3o(^) 

ao(r) = ao(fc)* , /3o(fc*) =/3o(A:)* (58) 



which can be easily solved to any order in k. Taking the 
chemical potential at half filling as the origin of energies 
and vp = 1, the first few terms are 

ao(fc) = ao,2|A:p + ao,3(fc^ + fc*^) + ... 

/?o(fc) = r +/3o,2A:' + /?o,3|fc|'fc* + ... (59) 

where all the constants are real. The eigenstates are still 
given by ([5S)) if one replaces the phase of k by the phase 
oif}*o{k). 



V. ELASTIC STRAINS AS EFFECTIVE 
ELECTROMAGNETIC FIELDS 

Here we show that the couplings of long wave- 
length phonons to electrons have some similarities — and 
differences — with those of the scalar and vector poten- 
tials for the electromagnetic field. We show, in particu- 
lar, that both in-plane and out-of-plane strains can mimic 
some of the effects of electric and magnetic fields. 

The minimal coupling prescription ki ^ ki + Ai on ([2|) 
gives 



H-ik) = 



$ -k- A 
-k* -A* $ 



(60) 



where A = A^ + iAy is the vector potential in complex 
notation and $ is the scalar (electric) potential. For the 
purposes of this section it is convenient to change from 
the circularly polarized modes of table I to linear modes 

Ui^—{Ux+lUy) , U2 = --J=[Ux - lUy) ^ 
Ql^ ^{Qx-lQv) , Q2^ ^{Qx+iQy) (61) 



Then, the acoustic hamiltonian (|49[) can be rewritten by 
using 

iV2{qQi + q*Q2) = '2.i{qxQx + lyQy) 
iV2qQ2 = 2i{qxQx - qyQy) - 2{qxQy + qyQx) 
iV2q*Qi = 2i{qxQx - QyQy) + 2{qxQy + qyQx) (62) 



This agrees with the hamiltonian given in reffi^ where the 
off-diagonal elements (second and third lines) were ob- 
tained by expanding a phonon-modulated hopping, while 
the magnitude of the diagonal element, known as 'defor- 
mation potential', was estimated in a nearly free electron 
model. Our derivation shows that both diagonal and off- 
diagonal terms are uniquely determined by symmetry. 

Taking Fourier transforms with iqi — > di and compar- 
ing with (|60p finally yields 



^Ki) = $(i^2) = 2ai{d,Q, + dyQy) 
A^Ki) = ~A4K2) = 2l3i{d.Qx - DyQy) 
Ay{Ki) = -Ay{K2) = -2l3iid^Qy + dyQ^) (63) 



where we have used (|35|) to obtain the electron-phonon 
hamiltonian around K2- These effective fields can be 
written in terms of elastic strains with the usual defini- 
tions 



dxQx 



dyQy , 2uxy = dxQy + dyQx (64) 



Thus we have the remarkable result that, around each 
Fermi point, static strains can mimic the effects of exter- 
nal electric and magnetic fields. 

We also see that the couplings of electrons near the 
two Fermi points are identical for the scalar potential, 
but differ by a sign for the vector potential. Thus, for 
static strains they will experience the same 'electrostatic' 
field, but opposite 'magnetic' fields given by 

B^iKi) = -B,iK2) = d,Ay{Ki) ~ dyA,{Ki) 
= -2/3i [{dl - dl)Qy + 2d,dyQ.,\ 

= -2fi\\ldxUxy + dyiUxX - Uyy)\ (65) 

We now address the following question: Is the (par- 
tial) identification between elastic strains and effective 
electromagnetic fields peculiar to graphene alone, or can 
we expect to find it in other systems? The first observa- 
tion is straightforward: The minimal coupling prescrip- 
tion k ^ k-\- A gives rise to the required type of interac- 
tions for the vector potential only if the electron system 
satisfies the Dirac equation. An ordinary non-relativistic 
equation would produce terms of the type k ■ A, where 
k is the electron momentum, unlike the leading phonon 
couphngs ((49|) . which are fc-independent. 

But even with Dirac points, further conditions have to 
be satisfied. The effective fields associated to acous- 
tic phonons are of the form diQj, where both di and Qi 
belong to the vector representation V . The vector poten- 
tial also belongs to V. In the continuum, V x V and V 
have opposite parities, and diQj and A can not transform 
equivalently. Thus, the possibility of describing elastic 
strains as effective fields is a lattice effect. More con- 
cretely, we need a lattice without inversion symmetry. 

But even on a lattice this identification can be partial 
at best. The reason is that the vector potential is odd un- 
der time reversal, whereas phonons and their derivatives 
are even. Since time reversal takes K to —K, as long as 
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—K is not equivalent to we can still hope for a partial 
identification valid around individual Fermi points. This 
will be possible only if the representations V and V x V 
of the little group have at least one irreducible rep- 
resentation in common, i.e., if V'^ contains the identity 
(or trivial) representation. 

We now discuss how these conditions are met by mono- 
laye graphene. The littele group is G^j = C^y with a 



vector representation that decomposes according to^ 



V ^ Ai{z)+E{x,y) 



(66) 



For in- plane modes only E is relevant, and the use of 
elementary group-theory techniquesiii^ yields 

ExE = Ai{d^Q^ + dyQy) + A2id^Qy - dyQ^) 

+ E{dyQy-d,Q,,d^Qy+dyQ,) (67) 

Here we recognize the L.H.S. of as the basis for the 
irreducible representations Ai and E. But the fact that 
the vector potential (Ax,Ay) and certain components of 
the strain tensor transform equivalently under the little 
point group Cs^ is not sufficient to guarantee that elas- 
tic strains can mimic a magnetic field around each Fermi 
point, for they could still couple to the two atoms in the 
unit cell with different signs. By (HU]), the two atoms are 
exchaged under G2T. Now, the vector potential is odd 
under G2 and T, whereas the the in-plane components of 
the elastic strain are even under both symmetries, and 
the two minus signs cancell each other. However, elec- 
trons at the two Fermi points see effective vector poten- 
tials which differ by the sign. On the other hand, the 
scalar potential $ belongs to the trivial representation 
Ai and is even under time reversal. As a consequence, 
the effective scalar potential takes the same sign on Ki 
and K2 = —Ki- 

This analysis can be extended to include the effects 
of so-called ripples or long wavelength deformations per- 
pendicular to the grapehene sheet. Such ripples make the 
two dimensional graphene sheet thcrmodinamically sta- 
ble and have been recently observed in individually sus- 
pended sheets?^. Although we have restricted ourselves 
to in-plane modes, acoustic out-of-plane modes can be 
easily incorporated in the long wavelength limit (see Sec- 
tion 6 for the possibility of a more general analysis): 

According to ([M]) a vertical displacement Qz belongs 
to the IR Ai of Gsi,. On the other hand, Qz is odd under 
reflection by the horizontal plane ah, whereas electronic 
wavefunction bilinears are necessarily even. This forces 
us to consider quadratic functions of Qz- By translation 
invariance along the z axis only derivatives of Qz are 
acceptable, and we have to consider quadratic functions 
of dxQz and dyQz- The basis for the IRs Ai and E are 
now given by 

Ai {{d^Qzf + [dyQzf) 
E {{dyQzf - {dxQzf,2dxQzdyQz) (68) 

The discussion after Eq. (|67p applies also in this case, im- 
plying that the couplings of electrons near the two Fermi 



points are identical for the scalar potential, but differ 
by a sign for the vector potential. Thus, the associated 
effective fields are given by 

$(i^l) = $(i^2) = 7 {{dxQzf + {dyQ.f) 
Ax{K^) = ~Ax{K2) - S {{dyQzf - (d^Qz)^) 

AyiKi) = -Ay{K2) - 25dxQzdyQz (69) 

where 7 and 5 are real coupling constants. 

We close this section by noting that frozen optical 
modes can also give rise to effective magnetic fieldsi^ 
Indeed, the substitution of ((6T|) in ((52|) gives 



(70) 



H{k,q,Q') = *x/2/3i (^ _^, ^Ji^ 

Q' + 



Q'y~^Q'x 



Comparing with the minimally coupled hamiltonians (|60p 
yields the identifications 



A',{K,) = -A',{K2) = P[Q'y 

4(ifi) = -A'y{K2) = -/3iQ^ (71) 

and the magnetic fields 

B'ziK^) - -B'ziK^) = -P'Ad^Q', + dyQ'y) (72) 

VI. DISCUSSION 

In this paper we have exploited the fact that, for many 
physical processes of interest in monolayer graphene, the 
wavevectors of electrons and phonons lie near points of 
high symmetry in the Brillouin zone. Even though the 
little groups for non-vanishing k and q are generically 
trivial, the proximity to points of high symmetry impose 
strong constraints on the series expansions of observables 
around them. This is analogous to the situation in Lan- 
dau's theory^ of second order phase transitions, where 
the dependence on the order parameter of observables in 
the low symmetry phase is determined by the space group 
of the high symmetry phase. Here, q and k play the role of 
order parameters, 'breaking' the symmetries of the F and 
Fermi points respectively. As a consequence, the method 
presented in this paper may be useful in other systems 
with Fermi points, such as multilayer graphene j^i^i^liS^ 
Semimetals with small electron and hole pockets around 
high symmetry points are also good candidates. 

Note that, in spite of the fact that we use tight-binding 
Bloch functions as our starting point, our results are 
model-independent. The reason is that only the sym- 
metry properties of the wavefunctions are used, and the 
general electronic hamiltonian ([57]) describes any doublet 
of states transforming according to the small representa- 
tion E of the little group. 

The relatively involved constraints (|40l4ip on the 
electron-phonon hamiltonian are reduced to extremely 
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simple conditions on the few functions which appear in 
the general solution in Section 4 — so simple, indeed, 
that they can be solved explicitly. In the process the gen- 
eral functions a and f3 are written as linear combinations 
of the symmetry- adapted coordinates. By dropping the 
linearity assumption we could extend our approach to 
multiphonon processes. One can easily check that the 
constraints (|40l4ip are still valid for nonlinear functions 
depending on several modes and momenta, with 

a{k,q,Q)^a{k,{qi},{Qi{qi)}) (73) 

and other obvious replacements. Note that at the non- 
linear level out-of-plane modes have to be included for 
consistency, and the whole symmmetry Dq^ rather than 
Cev has to be used. The analysis becomes more involved 
and will be the object of future work. 

Inspired by a mechanism first proposed in refj^ for 
the case of dislocations in multivalley conductors, effec- 
tive magnetic fields induced by elastic distorsions have 
been recently suggested as a way to explain the strong 
suppression of weak localization in graphenei^i^^ This 
differs from our results in Section 5 in that they consider 
effective gauge fields which are quadratic in the out-of- 
plane modes. Our analysis in the previous section shows 
that in-plane modes give rise to efffective fields which 
are linear in the strains. For physical applications of the 
treatment in Section 5 to dislocations and other deffects 
see TeiM 

Our method can be easily extended to accomodate ex- 
ternal fields. For instance, in the presence of an elec- 
tric field along the z axis, out-of-plane phonons will cou- 
ple linearly to the electrons. This may be an exter- 
nally applied electric field, or the effect of a substrate 
which breaks the symmetry under reflections on the sam- 
ple plane. The recently observed ripples in suspended 
graphene sheets^^ also break the reflection symmetry ah 
and can play the role of a 'background field', inducing 
new electron-phonon interactions. 
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APPENDIX A: TRANSFORMATIONS OF 
MODES AND WAVEFUNCTIONS 

The action of the generators of Cgt, on the electron and 
phonon wavevectors is given by the vector representation. 
This is simpler in complex notation 

C^q = Loq , aviq = q* , C2q = -q , OJ = e^"/^ 

(Al) 



where q = qx + i<ly 

Careful inspection of Table [1] shows that the generator 
acts identically on primed and unprimed symmetry- 
adapted modes, namely 

: ui{q) uj*ui{ujq) , U2iq) uJU2{Loq) (A2) 

where to = exp(27ri/3). The actions of CTui and C2, how- 
ever, are different for the two types of modes 

0-^,1 : Ui{q)^U2{q*) , U2{q) ^ Ui{q*) 

u[{q) ^ -u'2iq*) , u'.iq) ^ -u[iq*) 
C2 ■■ u^{q) ^ -u^{-q) , u[{q) u[{-q) (A3) 

Gorrespondingly, the symmetry-adapted coordinates Qi 
and Q'l in ([3]) have transformation properties 

gQ{gq)^T{g)Q{q) , 9Q\9q)=T\g)Q'{q) (A4) 

where T'(C2) = —T{C2) = 1 is the unit matrix, and 

T{Ct)^T'{Ct)^n^ °) ^^^^ 

Note that at the F-point the unprimed (primed) modes 
transform according to the two dimensional representa- 
tion El {E2) of Cei,. 

The action of the generators of the little group Cj,^ on 
the electronic wavefunctions is given by 

Cg^ : ^i+{k) cj$i+(wfc) , $2+(fc) w*$2+(wfc) 
^i-{k) — > tj*(f>i_(wfc) , ^2-{k) CJ$2-(t^fc) 
a^i : $i+(fc) ^ ^2+ik*) , $i_(fc) ^ <p2-{k*) (A6) 

The binary axis C2 connects the two Fermi points 

C2 : $i+(fc) ^ <^2-{-k) , $2+(fc) ^ 3>i-(-fc) (A7) 

Equations (jA6IA7p imply that, at the Fermi points 
{k — 0), the wavefunctions {<I'i±(0)} form the basis for 
the 4-dimensional IR E of the space group,— with star 

{Ku~Ki}. 

APPENDIX B: SYMMETRIES OF THE 
DYNAMICAL MATRIX 

Since according to (jl3p invariance under the generators 
of Cs^ impose identical constraints on D and D' , we write 
these only for D 

C+ : D{Luq) = nD{q)n* , : Diq*) = (T,i?(g)a, 

(Bl) 

The nondiagonal block M satisfies the same constraint 
under , but gets an additional minus sign under (jyi 

a„i : M{q*) = -axM{q)ax (B2) 
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Invariance under C2 implies that D and D' are even func- 
tions of q, whereas M is odd 

C2 : D{-q) = Dig) , D'{-q) = D' {q) , M(~q) = -Miq) 

(B3) 

Finally, time-reversal invariance ^ imposes the addi- 
tional constraint 

D{-q)^a,D*{q)a, (B4) 

which is also satisfied by D' and M. Combining this 
equation with (jBl|l - (IB3P gives 

D{q*) = D*{q) , D'{q*) ^ D'* (q) , M{q*) = M* [q) 

(B5) 

These constraints are solved by 

- ( nl) ) 

where a{q) and h{q) are even functions of q that satisfy 

a{ujq) = a{q) , biujq) = ojh{q) 
a{q*)=a{q) , h{q*)^b*{q) (B7) 

with identical solution for D'{q). For the off-diagonal 
block the solution is 

^'<') = ( 4'-) ) 

where and e{q) are odd functions of q that satisfy 
the following constraints 

d{ujq) — d{q) , e{ujq) = Lue{q) 
diq*)^d*{q) , e{q*)^e*{q) (B9) 

These imply the expansions 

d{q) = d3g'+rfV + 0(q^) 

e{q) - eiq + e3kl'g + 0(<z') (BIO) 

where all the constants are real. 



APPENDIX C: LITTLE GROUP CONSTRAINTS 
ON THE ELECTRON-PHONON HAMILTONIAN 

Substitution of p8|) into (j4T|) yields the following con- 
straints on the complex functions a and (3 

: a(wfc,wg,C3^Q) = a(fc,g, Q) 
P{u:k,u:q,C+Ct)=u*P{k,q,Ci) 
a„i : a(fc*,g*,cr„iQ) = C2Q) 

/3(r,(7*,a„iQ)=/?*(fc,g,C2Q) (CI) 



where the argument of Q is always the second one in 
the function, i.e., a{uk,Luq,C^Ql) actually stands for 
a{Lok,ujq,C^Cl{ujq)). The general solution to the con- 
straints can be found by writing 

a(fc,g, Q) = a{k,q,Q) + a' {k,q,Q') 
(3{k,q,Q) = (3{k,q,Q)+p'{k,q,Q') (C2) 

with 

a(fc,9,Q) = t{fl{k,q)Qiiq) + f2ik,q)Q2{q)) 
Pik,q,Q) = l{giik,q)Q,iq)+g2{k,q)Q2iq)) (C3) 

and equivalent expressions for a' and /?' in terms of 
primed variables. Then, substitution of a into the her- 
miticity conditions (|40l) gives 

/2(fc,g)--/i(fc,-g)* (C4) 

which implies 

a{k, q, Q) = z{f{k, q)QM ~ fik, -q)Q2{q)) (C5) 

with identical results for a'. Imposing the conditions (140]) 
on (3 and /?' gives 

9i{k,-q) ^ -gi{k,q) , g[{k,~q) ^ g[{k,q) (C6) 



Substitution into the little group constraints ()C1|) 
shows that the functions /, /', (72, and 52 must satisfy 
identical conditions, which we write only for / 

f{uk,uq)^ujf{k,q) , f{k*,q*)^ f{k,qy (C7) 

whereas the conditions on gi and g'l are 

gi{ujk,ujq) = gi{k,q) , gi{k*,q*) ^ gi{k,q)* (C8) 

The general solution for the electron-phonon hamiltonian 
is then given by (|42l43p . 
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